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^ ' The critical current modulation and the spontaneous flux of the vortex 

I states in corner Josephson junctions between Sr2Ru04 and a conventional 

s-wave superconductor are calculated as a function of the crystal orientation, 
^ I and the magnetic fleld. For Sr2Ru04 we assume two nodeless p-wave pairing 

states. Also we use the nodal /-wave states Big x and i?2g x E^, and 
^ ^ one special p-wave state having line nodes. It is seen that the critical current 
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I. INTRODUCTION 



During recent years, a renewed interest has been realized for unconventional superconduc- 
tivity in Sr2Ru04 (SRO) This material has a layered structure and is a two dimensional 
Fermi liquid, with similarities to ^He, that shows p-wave superfluidity. Knight-shift mea- 
surements suggest that unconventional superconductivity with pairing in the triplet channel 
is realized |^. The time-reversal symmetry is violated for this state. This is consistent with 
the discovery of internal magnetic flelds in the superconducting phase by yuSR measurement 
0. Also speciflc-heat measurements support the scenario of line nodes within the supercon- 
ducting gap as in the high- Tc cuprate superconductors |Q. Moreover thermal-conductivity 
measurements indicate the presence of horizontal lines of nodes The argument for un- 
conventional superconductivity in SRO is further supported by the variation of with 
nonmagnetic impurities 

The tunneling effect in ferromagnet-triplet-superconductor junctions can be used to de- 
termine the nodal structure of the pair potential and to distinguish the proposed pairing 
states. When the interface is normal to the ab plane, for two-dimensional order parameters, 
tunneling peaks are formed due to bound states that depend both on the surface orienta- 
tion and the quasiparticle trajectory angle |]^. For three dimensional order parameters the 
tunneling spectra is sensitive to the orientation of the interface |Q. 

The Josephson effect has been observed experimentally in bilayer junctions 
Pb/Sr2Ru04/Pb, where a sharp drop of the critical current Ic is observed below the tran- 
sition temperature of Sr2Ru04 |^. This anomalous effect originates from the competition 
between the positive Josephson coupling strength between the two singlet superconductors 
and the negative one due to the singlet-/triplet-superconductor coupling [0. Also in Ref. 



11 1 it is argued that the Josephson coupling between a conventional s-wave superconductor 



In and Sr2Ru04 is allowed in the in-plane direction, but not in the c axis. The conclusion 
is that the pairing state in that case has to be a nodeless p-wave state. 

In order to distinguish the pairing states with a two-dimensional order parameter, in 
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the present paper we study the static properties of a two-dimensional corner junction be- 
tween SRO and a conventional s-wave superconductor. This is done by mapping the two 
dimensional (2D) junction into two single junctions, connected in parallel, by introducing 
an extra relative phase in one part of this junction. The spontaneous flux and the critical 
current modulation of the vortex states with the junction orientation are calculated by solv- 
ing numerically the sine-Gordon equation. For SRO we shall assume three possible pairing 
states of two-dimensional order parameters, having line nodes within the Ru02 plane, that 
break the time-reversal symmetry. The first two are the 2D, /-wave states proposed by 



Hasegawa et ai, |jT2| having Big x and x Eu symmetry, respectively. The other one is 
called the nodal p-wave state and has been proposed by Dahm et al. [jl3[ , where the pairing 
potential has the form d = Aoz[sin(A;^a) -|- ism{kyay\. This pairing symmetry has nodes as 
in the x E^ case. Also we will consider two nodeless pairing states. One is the isotropic 
p-wave state and the other is the nodeless p-wave state initially proposed by Miyake and 
Narikiyo [|I^, both breaking the time-reversal symmetry. In addition we examine briefly 



several three-dimensional order parameters with horizontal lines of nodes. 

For each pairing state we derive simple arguments that could be used to characterize 
it. For example the junction orientation where the a, 6 axes point towards the junction 
interface will give finite critical currents only for the isotropic p-wave. Big x Eu states and 
zeros for the other states. Furthermore the suppression of the critical currents of the nodal 
state Big X Eu near the nodes can be used to distinguish it from the nodeless isotropic 
]9-wave state. Also the modulation of the Josephson critical current with the polar angle 
/?, within the ac plane, will probe the horizontal nodes for hypothetical three-dimensional 
order-parameter symmetries. The rest of the paper is organized as follows. In Sec. II we 
discuss the Josephson effect between a superconductor breaking time-reversal symmetry and 
a conventional s-wave superconductor. In Sec. Ill the geometry of the corner junction is 
discussed. In Sec. IV we present the magnetic flux for the spontaneous vortex states. In 
Sec. V the modulation of the magnetic flux and the critical currents with the orientation are 
presented. In Sec. VI several three-dimensional order parameters are considered. In Sec. 



VII the modulation of the vortex states with the magnetic field is discussed. In Sec. VIII a 
connection with the experiment is made. Finally, a summary and discussion are presented 
in the last section. 



II. JOSEPHSON EFFECT BETWEEN SRO AND A CONVENTIONAL 5- WAVE 

SUPERCONDUCTOR 

We consider the junction shown in Fig. |l|(a), between a superconductor A with a two- 
component order parameter n^, and a superconductor B with order parameter hb- The two 
superconductors are separated by an intermediate layer. The bulk order parameters ua and 
hb, near the interface, can be written as 

r^A = < (1) 

and 

nB = nBe"f'^, (2) 

where (I)a^, (t^Ay are the phases of the components of the order parameter ua, and </)b is the 
phase of the superconductor B. ua^, UAy are the magnitudes of the components of the order 
parameter ua- Then the supercurrent density can be written as 

J = Jba^ sin(0B - 0aJ + JsAy sin(0B - 0a J, (3) 

where 

Jba^ ~ UA^riB, 

(4) 

JsAy ~ riAyTlB- 

We have to comment here that the coupling between a triplet and a singlet supercon- 
ductor is forbidden by the orthogonality of the spin part of their wave functions. However, 
in our situation there exists a finite overlap due to spin-orbit coupling |T3. In the presence 



of spin-orbit coupling the Cooper pairs of different parities will be mixed at the interface 



between the singlet and the triplet superconductor resulting in a direct Josephson coupling 
between them [0]. We define (p = (ps — (pA^ as the relative phase difference between two 
superconductors. We consider the case where the intrinsic phase difference within supercon- 
ductor A is (f)Ay — <pA^ = tt/2. In this case the order parameter is complex and breaks the 
time-reversal symmetry. Then the supercurrent density can be written as: [ p!7| ] 

J(0) = Jesin(0 + 0J, (5) 

with 



tan-i -4^ Jba. > 0, 



JbAi 



n + tan-i ^ Jsa. < 0, 



(6) 



and 



Jc - y^BAy + Jba^- (7) 

We consider the following cases for the pairing state of SRO having line nodes. 

(a) In the first 2D /-wave state Big x Eu-, the magnitudes of the order parameters ha^ 
and riAy in Eq. (^) are defined as ua^ = cos(2a) cos(q;), riAy = cos(2q;) sin(a). In these 
formulas a denotes the angle between the interface and the a axis of the crystal. This state 
has nodes at the same points as in the dx2_y2-wdNe case. 

(b) For the second 2D /-wave state i?2g x E^-, tla^ = sin(2a) cos(a;), hAy = sin(2a) sin(a). 
This state has nodes at 0, 7r/2, tt, and 37r/2, and has also been studied by Graf and Balatsky 



(c) This is the case of a nodal p-wave superconductor, ua^ = 1/sm sin[7rcos(a)], nAy 



1/sm sin[7r sin(a)]. We use here the same normalization proposed by Dahm et al, fTS 
sm = sin(7r/-\/2) = 1.125, where the Fermi wave vector is chosen as kpa = vr, in order 
to have a node of the order parameter in the Fermi surface. This state has nodes as in the 
X Eu state. The corresponding nodeless form was initially proposed by Miyake, and 
Narikiyo and will be considered as a separate case. Also we will examine the following 



two pairing symmetries, which are nodeless 



d) In the case of a nodeless p-wave superconductor, proposed by Miyake and Narikiyo 
= 1/sa/ sin[i?7r cos(a;)], fiAy = I/sa/ sin[i?7rsin(Q;)], Sm = sm{TC / y^) = 1.125, 
and R = 0.9. This state does not have nodes. 

e) In the isotropic p-wave case, = cos(q;), UAy = sin(a). This pairing state does not 
have nodes either. 



We consider the corner junction shown in Fig. |T|(b) between a triplet superconductor 
breaking the time-reversal symmetry and an s-wave superconductor. In the following the 
Josephson coupling is restricted in the ab plane since we assume 2D order parameters. If the 
angle of the a axis with the interface in the x direction is a, then the corresponding angle 
in the y direction will be 7r/2 — a. Using this constraint we map the two segments of this 
corner junction, each of length L/2, into a one-dimensional axis shown in Fig. |l|(c). In this 
case the two-dimensional junction can be considered as being made of two one-dimensional 
junctions described in Sec. II connected in parallel. The characteristic phases of the two 
parts of the one- dimensional junction 0ci and 0c2 depend upon the angle a. We call this 
junction frustrated since the two segments have different characteristic phases (pci and 0c2- 
The fabrication details of corner junctions or SQUID's, between sample faces at different 
angles, can be found in Refs. [|l^ and [pO[] . 

The superconducting phase difference (p across the junction is then the solution of the 
sine- Gordon equation 



with the boundary conditions 



where (pdx) = 0ci(0c2) in the left (right) part of the junction, and I"" is the bias current in 



III. THE JUNCTION GEOMETRY 




(8) 
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the overlap geometry. The length x is scaled in units of the Josephson penetration depth 
given by 



A. 



V SnedJco' 

where d is the sum of the s-wave, and mixed wave, in-plane penetration depths plus the 
thickness of the insulator layer. The different solutions obtained from Eq. (|^) are classified 
with their magnetic-flux content 

<^ = ^{<Pr-<Pl), (10) 

where 4>r(l) is the value of (j) at the right (left) edge of the junction, in units of the flux 
quantum $o = hc/2e. 

To check the stability of solutions we consider small perturbations u{x,t) = u{x)e'^* on 
the static solution 0(x), and linearize the time-dependent sine-Gordon (s-G) equation to 
obtain 

dPvix^ ~ 

\- JcCOs[0(x) + (t)c{,x)]v{x) = \v{x), (11) 

where A = — s^, with the boundary conditions 

dv , 



dx 



\x=0,L 



0. (12) 



It is seen that if the eigenvalue equation has a negative eigenvalue the static solution (f){x) 
is unstable. 



IV. SPONTANEOUS VORTEX STATES 

First let us consider the case where the two one-dimensional junctions between a triplet 
superconductor breaking the time-reversal symmetry and an s-wave superconductor, each of 
length L/2, described in Sec. II, are uncoupled. Then for < x < L/2 the stable solutions 
for the s-G equation are (f){x) = — 0ci + 2ni7r, where rii = 0, ±1, ±2, while for L/2 < x < L 
the stable solutions for the s-G equation are 0(x) = — 0c2 + 2n2vr, where n2 = 0, ±1, ±2, 
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and 0ci, 0c2 are the relative phases in each part of the junction due to different orientations. 
When the frustrated junction is formed, and we consider the above junctions in parallel, the 
phase (j) is forced to change around x = L/2 to connect these stable solutions. This variation 
of the phase along the junction describes the Josephson vortices. The flirx content of these 
states (in units of $o) is [0 

$ = - 0(O)]/27r = (-0,2 + 0ci + 2n7r)/27r, (13) 

where the n value [n = ni — n2 = 0, ±1, ±2, ...) distinguishes between solutions with different 
flux content. We will concentrate on solutions called modes with minimum flux content, i.e., 
n = 0, 1, — 1. Their magnetic flux in terms of (pci, 4>c2 is shown in Table |. Generally the flux 
content is fractional, i.e., is neither integer nor half-integer, as a consequence of the broken 
time-reversal symmetry of the problem. 

V. SPONTANEOUS MAGNETIC FLUX AND CRITICAL CURRENTS 

In this section we will consider the magnetic flux and critical current modulation with the 
azimuthal angle a, for the various two-dimensional pairing states for the SRO. In general the 
flux quantization and the critical currents are influenced by the nodal structure of the order 
parameter. The critical current 1°^ which depends mainly on the magnitude of the order 
parameter vanishes at the nodes of the order parameter. The magnetic flux is determined 
from 0ci;0c2; ciud has almost the same structure in the isotropic p-wave. Big x and 
X Eu pairing states, while it is different in the nodal p-wave and the nodeless p-wave 
cases. 

A. Isotropic p wave 

For the isotropic p-wave case we present in Fig. |^(a) the variation of 0ci,0c2, with 
the interface orientation angle a. The magnetic flux can be calculated analytically for the 
spontaneous vortex states 0, —1, 1. The corresponding expressions are seen in Table |I|. The 
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spontaneous magnetic flux as a function of the angle a is presented in Fig. |](b). At a = 
we have the existence of three possible spontaneous vortex states. As we change the angle 
a, the magnetic flux changes, keeping the distance between the modes at a constant value 
equal to a single flux quantum. The modes for a = are displaced to quadratic values [i.e., 
(1/4 + n)$o] of the magnetic flux in contrast to the s-wave case where the flux is quantized 



in integer values 0, —1, 1 P2|l . 



In particular the vortex solution in the n = mode (solid line) contains 1/4 of a fluxon 
for a = 0, and as we increase the angle a towards 7r/4 it continuously reduces its flux, i.e., 
it becomes flat exactly at a = 7r/4 and then it reverses its sign and becomes an antivortex 
with exactly opposite flux content at a = it/2 from that at a = 0. In addition we have 
plotted in Fig. |^(a) the phase distributions for the mode n = in different orientations 
a = 0,7r/4,7r/2. The transition from the vortex to the antivortex mode as the orientation 
changes is clearly seen in this figure. Note that the solutions in this mode remain stable 
for all junction orientations. This is seen in Fig. ^(c) where we plot the lowest eigenvalue 
(Ai) of the linearized eigenvalue problem as a function of the angle a. We see that Ai > 0, 
denoting stability for all values of the angle a in this mode. 

Let as now examine the solution in the n = —1 mode [dashed line in Fig. 0(b)]. We see 
that at a = it has negative flux, equal to —0.75, and is stable. As we increase the angle 
a it decreases its flux to a full antifluxon when the orientation is 7r/4 and then to flux of 
—1.25 when a reaches 7r/2. As seen in Fig. |^(c) this solution becomes unstable at a = 7r/4 
where a full antifluxon enters the junction. The phase distribution at this point is seen in 
Fig. 1(b) (dotted line). 

Finally the solution in the n = 1 mode [long-dashed line in Fig. 0(b)] contains 1.25 in 
flux at a = and is clearly unstable. It becomes stable at a = 7r/4 [see Fig. ^(c) (long- 
dashed line)], where a full fluxon enters the junction, as seen in Fig. 0(c) (dotted line). At 
a = 7r/2 it contains 0.75 in flux. 

In Fig. 0(d) we plot the overlap critical current per unit length 1°"" as a function of the 
angle a, at H = 0, for the n = 0, —1, 1-mode solutions. The Josephson critical current 
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density is kept constant with the orientation Jc = 1. Let us consider the situation where 
the junction contains a solution in the mode n = 0, at a = 0, when the net current is 
maximum. The spatial variation of is described by a fractional vortex displaced from the 
corresponding distribution at zero current which is around 7r/4 [see Fig. 3(a)] . The current- 
density distribution as seen in Fig. ^(a) (solid line) corresponds to a positive fluxon, and at 
the maximum current is flat, close to unity, and has a small variation around the junction 
center giving rise to the large value of the net current, seen in Fig. 0(d). At a = tt/A 
the flat phase distribution corresponding to the n = solution seen in Fig. §(a) at zero 
current is displaced towards the value Sir /A at the maximum current. The current-density 
distribution as seen in Fig. 1(a) (dotted line) is flat, giving a net current equal to 1 which 
is the maximum value over the range < a < n /2. 

For the —1 solution at a = the phase represents an antifiuxon which is being displaced 
at the maximum current from the corresponding distribution at zero current seen in Fig. 
^(b). The corresponding current density represents an antifiuxon, and has a characteristic 
variation along the junction, taking positive and negative values, as seen in Fig. ^(b) (solid 
line), so that the net current is suppressed compared to the mode. At a = 7r/4 just one 
antifiuxon enters the junction and the current distribution is symmetric around x = L/2, as 
seen in Fig. §(b) (dotted line), leading to zero total current at this angle. 

The critical current 1°^ for the mode is finite for each orientation a due to the nodeless 
form of the order parameter. This is in marked contrast to what we see in the singlet pairing 
states that describe the high-Tc cuprates where the critical current 1°"" for the different modes 
is suppressed near a = tt/4 following a similar variation to the Josephson critical current 
density Jc(ct) Also the fractional values of the magnetic fiux are originated from the 
two-component order parameter that breaks the time-reversal symmetry. 
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B. Big X 



In the Big x i?„-wave pairing state the (pci, 4>c2 and also the magnetic flux seen in Fig. |^(a) 
and Fig. ^(b), respectively, have the same form as in the isotropic p-wave case. This happens 
because the relative phase 0c in Eq- ® depends on the ratio of the real and imaginary 
parts of the order-parameter magnitude in superconductor SRO. The Big component cos2q; 
multiplies both parts and is canceled out. On the contrary only the Big part of the order 
parameter contributes to the Josephson critical current density Jc, which varies as | cos(2q;)| 
as seen in Fig. ^(a) (dashed line), while the part is canceled. It has nodes at a = ±7r/4. 
For this reason we choose the interval —it /A < a < 7r/4 to study the Big x Eu pairing state. 

For a = the mode n = contains negative magnetic flux of —0.25 as seen in Fig. |^(b), 
the mode 1 contains positive magnetic flux equal to 0.75 and is stable, while the mode —1 
has —1.25 in flux and is unstable. The variation of the phase of the vortex solutions at zero 
current and the corresponding current density at maximum current for a = are seen in 
Fig. ^. Differently from the isotropic p-wave case, the mode with negative magnetic flux 
corresponds to the highest critical current. 

Close to ±7r/4 the variation of the magnetic flux for the modes 0, —1, 1 is influenced from 
the suppression of the Jc which goes to zero at ±7r/4 and the magnetic flux acquires some 
curvature at these points as seen in Fig. ^(b). Due to this abrupt variation the numerical 
procedure to follow modes —1,1 stops at the point where the flux changes more rapidly. 
This is the reason why the magnetic flux for the modes —1,1 does not reach the nodes of 
the Jc at a = ±7r/4. 

Also the critical currents are influenced by the Jc- The critical current for the mode 
vanishes at a = ±7r/4 where the order parameter has nodes. Also the maximum of the 
critical current tends to move towards a = where Jc has its maximum value. On the other 
hand the phase prefers the orientation a = —tt/A where it is flat. These two competitive 
factors shift the maximum current for n = to the left of a = 0. The critical currents 
for the other two modes are suppressed close to a = ±7r/4 where the order parameter has 



11 



nodes. In particular for mode 1 the magnetic flux is smaller compared to mode —1 and as 
a consequence the critical current is higher. In the interval we study, mode —1 is unstable 
while modes 1,0 are stable as seen in Fig. ^(c), where the lowest eigenvalue Ai versus a is 
plotted. 

C. B2g X Eu 

In the x case the 0ci, 0c2 seen in Fig. 0(a) have the same form as in the isotropic p 
wave case. Here the i?2g factor of the order parameter sin 2a is canceled out in the calculation 
of the relative phase 0c(a;). The i?2g part of the order parameter influences the Josephson 
critical current Jc, which varies as | sin(2a)| as seen in Fig. |^(a) (dashed line) and has nodes 
at 0, 7r/2. Due to this form of the Jc we have nodes at the critical current l'^ versus a at the 
points 0,71/2, as seen in Fig. |3(d). The modes 1, —1 have zero critical current at the angle 
a = 7r/4 where a full fiuxon or antifiuxon enters into the junction. In addition they become 
zero close to or 7r/2 due to the vanishing of at these angles. Note that close to or 
7r/2 the form of influences also the magnetic flux at zero current which acquires a small 
curvature at these points. The mode has its maximum close to a = in the B\g x Eu state 
while in x Eu the critical current at a = becomes zero. In the isotropic p-wave case 
the critical current is finite at a = since the order parameter for this state is nodeless. 

D. Nodal -p wave 

For the nodal p-wave case the 0ci, 0c2 are interchanged compared to the isotropic p-wave 
and B\g x E^, cases and also their linear form is lost. As a result the slope of the magnetic flux 
versus a seen in Fig. ^(b) has changed from negative in the isotropic p-wave and B^g x E^ 
cases to positive here. Now mode 1 is stable in the interval < a < 7r/4 and mode — 1 
is stable in the interval 7r/4 < a < it/ 2, contrary to the isotropic p-wave case as seen in 
Fig. ^(c). The Josephson critical current density Jc has the same nodal form as in the 
B2g X Eu case, with nodes at 0, 7r/2, as seen in Fig. ||(a), since the nodal form of the order 
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parameter is the same. The critical current for mode versus a follows the same form as 
in the x case and it is zero at 0, 7r/2 and maximum at a = 7r/4, where the phase is 
uniform. Also the variation of the flux versus a seen in Fig. |](b) is quite close to zero for 
mode 0, and as a consequence the 1°^ vs a for mode is more broadened compared to the 
corresponding mode in the -B23 x Eu case. The mode —1(1) becomes zero at a = 7r/4 where 
a full antifluxon (fluxon) enters the junction. For a < tt/4, mode 1 has less magnetic flux in 
absolute value than does mode —1 and its critical current 1°'" is higher. On the other hand 
for a > it/ 4, mode 1 has more magnetic flux in absolute value than does mode —1 and its 
critical current 1°"" is lower. Note that modes —1, 1 stop abruptly at the angle where the 
variation of the flux is more abrupt 

E. Nodeless p wave 

For the nodeless p-wave case the variation of 4'ci,4'c2 is more complex. It is similar to 
the nodal p-wave case except at values close to 0,7r/2 where it changes slope. As a result 
the magnetic flux versus a varies more closely to integer values and is almost flat as seen in 
Fig. ^(b). Also within the interval < a < 7r/2 modes —1, 1 change from stable to unstable 
for more than one time as seen in Fig. P(c) while mode is stable for all values of a. The 
Jc has a nodeless form as seen in Fig. |](a), i.e., it has no zero values. As a consequence the 
critical current for mode is finite over the angles a and has its maximum values at a = 7r/4 
where the phase is uniform. At this point a full fiuxon or antifluxon enters the junction and 
the critical current for modes —1,1 is zero. Also the magnetic flux for modes —1,1 varies 
close to values —1, 1 and as a consequence their critical currents are almost coincident. 

VI. 3D ORDER PARAMETERS WITH HORIZONTAL LINES OF NODES 

In the preceding sections we considered purely two-dimensional order parameters. In 
this section we consider the case where the order parameter is three dimensional. Recent 
measurement of thermal conductivity under a rotating in-plane magnetic fleld by Izawa et 
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al. has shown a fourfold pattern with a minimum when H \\ [110] and suggests that the 
nodes of the gap are on the hne with /c^=const as discussed in Ref. |12|. So the question 
arises also whether the results obtained in the present paper remain valid for such parings. 

Several three-dimensional order parameters have been proposed in the literature: 
Hasegawa et al. [O have proposed the pairing states with order parameter d = AQz{kx + 



iky) cos{ckz) and d = Aoz[sin(aA;^/2)+i sin(a/2)] cos(cA;2/2) withe being the lattice constant 
along the c axis. These states have horizontal lines of nodes at k^ = ±7r/2c and kz = ±vr/c, 
respectively, and break the time-reversal symmetry. Also Won and Maki have proposed 



an /-wave model with a superconducting order parameter d = Aoz{kx + iky^kz- This state 
has horizontal lines of nodes where the basal plane kz = intersects the Fermi surface. 

In all the above pairing states the Josephson current within the ab plane is canceled, since 
the order parameter changes sign in the kz direction symmetrically around zero. However it 
is possible to have finite Josephson current when the Fermi surface is not exactly cylindrical 
but is corrugated as observed by the angle dependence of the de Haas-van Alphen effect 
P5| . In that case the cancellation would not be complete. Another possibility to have finite 
Josephson current within the ab plane is to have a three-dimensional order parameter that 
is mixed with a two-dimensional one as suggested also in Ref. Then for the Josephson 
effect within the ab plane, discussed in preceding sections, the quasiparticles will experience 
an effective order parameter that is quasi-two-dimensional and the modulation of the critical 
current and magnetic flux with the misorientation angle within the ab plane will show the 
nodal profile of the order parameter within the ab plane discussed previously. 

In an alternative geometry, the modulation of the critical current and spontaneous flux 
with the polar misorientation angle within the ac plane, i.e., when the azimuthal angle 
is set to zero a = 0, would reveal the angular anisotropy of the order parameter along 
the c axis. In that case the effective order parameter felt by the quasiparticles is two 
dimensional. Moreover for that orientation the two-component order parameter is reduced 
to a one-component one which has a sign change but does not break time-reversal symmetry. 

For example we consider the three-dimensional order parameter d = Aoz(/c2,' + 
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iky) cos{ckz). For a spherical Fermi surface and a = 0, it is reduced to the two-dimensional 
order parameter d = Aqz sin/5cos(7r cos/3). We plot in Fig. |l3(a) the order parameter for 
adjacent crystal edges versus the polar angle (3. It is seen that the order parameter changes 
sign, in a restricted angle regime, at adjacent faces of the crystal. Also the characteristic 
phases 4'ci,4>c2 differ by tt as seen in Fig. ^n|(b). That would, for example, lead to half- 
integer spontaneous flux states. Note that the Josephson critical current densities as seen 



in Fig. |TO|(c) are not equal for adjacent junction edges. However the asymmetry in the 
current densities will not affect that spontaneous flux. The modulation of the current with 
the polar angle within the ac plane will be highly anisotropic for the three-dimensional order 
parameters. It will be suppressed for ck^ = tt/2 for the pairing state A oc coscfc^. 

Since all states considered here have d || z the Josephson current along the c axis between 
the singlet and triplet superconductors is not allowed In fact the spin-orbit coupling 
strength of the triplet superconductor would give rise to the orientational dependent Joseph- 
son coupling between the SRO and the conventional s-wave superconductor which completely 
vanishes along the c axis. However the Josephson coupling is possible when the d vector is 
not exactly parallel to the current. In that case we restrict to the corner junction in the ac 
plane where the c axis is tilted within the ac plane for certain polar angle. For the pair- 
ing state A oc kz the singlet-superconductor /triplet-superconductor/singlet-superconductor 
junction along a direction slightly different from the c axis can be used to probe the pairing 
symmetry from the sign change of the Josephson current at opposite faces of the junction 

Pi- 



VII. MAGNETIC FIELD 

We now examine the influence of the magnetic field on the spontaneous vortices, for the 
a6-plane Josephson effect, for fixed junction orientation a = 0, where most experiments on 
corner junctions have been performed. The critical current is finite, at a = for the pairing 
states isotropic p wave and Big x Eu, while it is zero or very suppressed for the other pairing 
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states. Hence we will restrict ourselves only to the isotropic p wave and Big x at a = in 
this section. The modes extend to a large interval of the magnetic field, which is almost the 
same. The instability at the extremum values of the magnetic field is due to the interaction 
of the fiux entered from the edges with the intrinsic fiux. The separation of the modes in 
magnetic fiux is kept constantly equal to $o as in the s-wave case, but is displaced compared 
to the s-wave case by an amount that is equal to the intrinsic fiux. 

We plot in Fig. |ll|(a) the overlap critical current 1°^ versus the magnetic field H, for a 
junction with length L = 10, for angle a = 0, and for the isotropic p-wave case. The mode 
0, at if = 0, contains positive fiux $ = 0.25 and has greater critical current than mode 
— 1, which at H = contains negative magnetic fiux and is equal to —0.75. So the mode 
with positive magnetic fiux corresponds to the highest critical current. The maximum I"'" 
for mode occurs to the right of H = 0. Note that the displacement of modes 0, —1, 1 in $ 
from their positions in the s-wave case is small and thus the 1°^ vs H diagram is similar to 
the s-wave case. 

In Figs. |T^(a) and |T2|(b), the overlap critical current versus the magnetic field H 



and the spontaneous fiux vs magnetic field is plotted for the Big x case. The mode 0, 
a.t H = 0, contains negative fiux $ = —0.25 and has greater critical current than mode 1, 
which at H = contains positive magnetic fiux and is equal to 0.75. So differently from 
the isotropic p-wave case, the mode with negative magnetic flux corresponds to the highest 
critical current. 

For the rest of the pairing states the junction orientation a = gives zero or relatively 
small values for the Josephson critical current density Jc(«), which means that the critical 
current the junction can carry will be zero or very suppressed compared to the previous cases. 
For the ac-plane Josephson effect for three-dimensional order parameters it is possible to 
have a sign change of the order parameter at adjacent faces of the corner junction that 
would, for example lead to an anomalous Fraunhofer pattern with "dip" at zero magnetic 
fleld. The asymmetry in the Josephson critical current densities will only cause the "dip" 
to be shallow and will maintain the symmetry of the pattern. 
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VIII. EXPERIMENTAL RELEVANCE 



For the a6-plane Josephson effect for two-dimensional order parameters, the measurement 
of the magnetic interference pattern at a = where the a and h crystallographic axes are 
perpendicular to the junction's faces will not give any finite value for the critical current 1°^ 
for the pairing states i?2g x Eu or the nodal p wave since the value a = is a node of the 
order parameter for these symmetries. For the nodeless p-wave case the Josephson critical 
current density is small at a = and so is the critical current. So the observation of finite 
critical current 1°^ in corner junctions at a = would be an indication for isotropic p-wave 
or Big X Eu pairing states where the value a = is not a node. Furthermore the suppression 
of the critical currents of the nodal state Big x E.^ near the nodes, i.e., at a = ±7r/4, can be 
used to distinguish it from the nodeless isotropic p-wave state, for which the critical current 
for no-flux mode does not become zero. 

On the other hand for the ac-plane tunneling, the orientation where ckz = 7r/2 will 
suppress the Josephson critical current for the pairing state A oc cos ck^ since it is a node 
for the order parameter. Also the singlet /triplet /singlet junction along the direction slightly 
different from the c axis will probe the pairing state A oc /c^ from the sign change of the 
Josephson current at opposite sides of the junction. 

Another experiment that we propose here to discriminate between the candidate's pairing 
states consists of a SQUID between a Sr2Ru04 thin film patterned into a circle with a series 
of Nb-Au-Sr2Ru04 edge junctions spaced at a fixed angle step. The tunneling directions can 
be defined lithographically and patterned by ion milling of a c-axis-oriented film to study the 
a6-plane Josephson effect or a 6-axis-oriented film to study the ac-plane Josephson effect. The 
measurement of the critical current versus the angle a, which mainly probes the magnitude 
of the order parameter, is expected to show large anisotropy with the angle a. An analogous 
experiment has already been performed in the SQUID geometry for YBa2Cu307 thin films 

A number of complicating factors are involved in the interpretation of the experiments 
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with corner junctions that could modify the anticipated results. These include the asymme- 
try in the Josephson critical current densities, the flux trapping, and the sample geometry. 
Also the study of the direction-dependent Josephson effect will be influenced by the rough- 
ness of the interface due to the existence of Ru lamellas, and also by the difficulty in cleaning 
and polishing a crystal at angles different than and 7r/2. However these effects are not 
sufficiently large enough to change the qualitative picture of the experiments. 

IX. CONCLUSIONS 

We followed the modulation of the vortices that are spontaneously formed in a corner 
junction between a triplet superconductor breaking time-reversal symmetry and an s-wave 
superconductor with the junction orientation angle. The nodal structure of the order param- 
eter is reflected in the Josephson critical current density, and influences the critical current 
versus a diagram and the magnetic flux of the spontaneous vortex states that are formed at 
the corner junction. The critical current vanishes at the angles where the order parameter 
has nodes. In addition it becomes zero at the angles where a full fluxon or antifluxon enters 
the junction. For the nodeless pairing state mode does not vanish, and has the nodal 
proflle of the corresponding order parameter. The presence of spontaneous magnetic flux 
can be used to distinguish between different symmetry states. 

The external magnetic field for the junction orientation a = where the a axis points 
towards the junction edge will only influence the states in which their order parameter is 
flnite at a = 0, i.e., isotropic p wave and B\g x E^. The critical current and the effect of the 
magnetic field will be zero or very suppressed for the other cases. For three-dimensional order 
parameters the Josephson critical current modulation, which mainly probes the magnitude 
of the order parameter with the polar angle, is expected to be highly anisotropic for the 
ac-plane Josephson effect. 

In corner junctions with broken time-reversal symmetry, the spontaneus flux and the 
magnetic interference pattern, even in the short junction limit, are influenced by the mag- 
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nitude of the order parameter and also the difference in the intrinsic phases 0ci and 0c2, at 
adjacent junction edges. This difference uniquely characterizes the pairing symmetry and 
can be used without any ambiguity to probe anisotropic pairing states with sign change of 
the order parameter between orthogonal directions in k space. 
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TABLE I. The magnetic flux (<I>) in terms of 0^1, (f)c2 for the spontaneous solutions that exist 
in a corner junction between a triplet superconductor with time-reversal broken symmetry and an 
s-wave superconductor {(f)ci,<f>c2 is the extra phase difference in the two edges of the corner junction 
due to the different orientations, of the a axis of the triplet superconductor). We present only the 
minimum flux states n = 0, — 1, 1. 

Vortex state Magnetic flux ($) 

(-</)c2 + <^ci)/27r 

1 {-4>c2 + -^ci + 27r)/27r 
-1 (-<^c2 + -^ci - 27r)/27r 



TABLE IL The magnetic flux ($) in terms of a for the spontaneous vortex states that exist in 
a corner junction between a triplet superconductor with isotropic p-wave symmetry and an s-wave 
superconductor. 

Vortex state Magnetic flux ($) 

0.25 - a/TT 

1 1.25 - a/TT 
-1 -0.75 -a/TT 
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FIG. 1. (a) A single Josephson junction between a triplet superconductor A with broken 
time-reversal pairing symmetry and B, with s-wave symmetry. The angle between the crystalline 
a axis of A and the junction interface is a. (b) The geometry of the corner junction between 
superconductors SRO and s wave, (c) The mapping of the corner junction geometry into the 
one-dimensional axis of length L. For < < L/2 the angle between the a axis and the interface 
is a while for L/2 < x < L the corresponding angle is 7r/2 — a. 
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FIG. 2. (a) The variation of (pci, 4'c2 with the angle a. For each value of a there exist three 
vortex solutions, (b) The spontaneous magnetic flux $ as a function of the angle a, for the 0, —1, 1 
mode solutions as solid, dashed, and long-dashed lines respectively, (c) The lowest eigenvalue Ai 
as a function of the angle a, for the 0, —1, 1 mode solutions, (d) Overlap critical current density 
1°^ versus the angle a for a junction of length L = 10. The pairing symmetry of the triplet 
superconductor is isotropic p wave. 
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FIG. 3. The phase distribution of the vortex solutions (a) n = 0, (b) n = —1, and (c) n = 1, 
at a = 0, 7r/4, 7r/2; for a corner junction of isotropic p-wave and s-wave superconductors, with 
length L = 10, and zero overlap external current 1°^ = 0. 
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FIG. 4. The current density distribution J(x) of the vortex solutions (a) n = 0, (b) n = — 1, 
and (c) n = 1, at a = 0, '7r/4, vr/2; for a corner junction of isotropic p-wave and s-wave supercon- 
ductors, with length L = 10, and maximum overlap external current 1°"" . 



26 




27 



-6.0 




FIG. 6. (a) The phase distribution of the vortex solutions at a = 0, for a junction with length 
L = 10, and zero external current I""" = 0. (b) The corresponding current density at the maximum 
external current. The pairing symmetry of the triplet superconductor is Big x Eu. 
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FIG. 10. (a) The effective two-dimensional order parameter d = zsin/3cos(7rcos/3) versus the 
polar angle /? for adjacent junction edges as solid and dashed lines, (b) The characteristic phase 
0ci (0c2) as a solid (dashed line), (c) The Josephson critical current densities Jd, Jc2 for adjacent 
junction edges as solid and dashed lines. 
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FIG. 11. (a) The variation of the overlap critical current with the magnetic field H for 
the three vortex solutions n = 0, 1, —1. (b) The spontaneous magnetic flux <I> as a function of the 
magnetic field H, for the 0, —1, 1 mode solutions, for a junction of length L = 10, and the angle 
a = 0. The pairing symmetry of the triplet superconductor is isotropic p wave. 
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